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Abstract 


Using the nonperturbative Schwinger-Dyson equation, we show that chiral 
symmetry is dynamically broken in QED at weak gauge couplings when an 
external magnetic field is present and that chiral symmetry is restored at 
temperatures above Tc — ■^\/2'K\eH\, where a is the fine structure constant 
and H is the magnetic field strength. 
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I. INTRODUCTION 


Two of us have recently proposed a method to study dynamical chiral symmetry breaking 
in gauge theories in the presence of an external held using the Schwinger-Dyson equation 
approach [|^. It is the purpose of this paper to describe the details of the methodology, 
using the case of a uniform magnetic held as an example. We also show how to adopt our 
formalism to study hnite-temperature ehects. 

We use quantum electrodynamics (QED) as our model gauge theory and consider chiral 
symmetry breaking in the presence of a constant external magnetic held. We introduce the 
formalism in Section II and derive the Schwinger-Dyson equation for the fermion self-energy 
in the quenched, ladder approximation, Eq. (|3^) . Using an approximation suitable for weak 


gauge couplings, we derive in Section III an approximate gap equation, Eq. (|5^) , from which 
the infrared dynamical fermion mass, Eq. , is obtained. Our result is consistent with that 

found by Gusynin, Miransky, and Shovkovy p[, who used a diherent approach. We show 
in the appendix how our formalism can be applied to the approach of Gusynin et ai, and 
establish the existence of the Nambu-Goldstone boson of the spontaneously broken chiral 
symmetry. 

It has been suggested in the literature [Q that the chiral symmetry breaking solution 
may hnd applications in the electroweak phase transition during the early evolution of the 
universe. To verify this possibility, it is necessary to take into account the thermal conditions 
present in the early universe. Our formalism makes it easier to study such hnite-temperature 
effects. This is discussed in Section IV where we obtain an estimate of the critical tempera¬ 
ture for chiral symmetry breaking, Eq. ([70|) . Our result indicates that, in order for the chiral 
symmetry breaking solution found here to be relevant for the electroweak phase transition, 
an unacceptably large magnetic held must be present at the time of the phase transition. 
We oher our conclusions in Section V. 
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II. FORMALISM 


Let us consider chiral symmetry breaking in QED in the presence of a static, external 
electromagnetic held. The Green’s function that describes the motion of a fermion with 
electric charge e in such an external held satishes the equation, 

7 ■ n(x)GA(x,|/) + J d‘^x'M{x,x')GA{x',y) = -y), (1) 

where n^(a;) = —id^ — eA^^^{x)^ and M{x,x') is the mass operator M in the coordinate 
representation: M{x,x') = {x\M\x'). As pointed out by Ritus |^, M is a scalar 7-matrix 
function of the fl^ and the and for constant 

M = {F,Mn",l 5 F^uFn- ( 2 ) 

In other words, for uniform external helds, only four independent 7-matrix valued scalars 
can be formed out of 11 ^ and F^^, as listed in Eq.(|D, where F^’^ = ^e^^'^'^Fxr- Furthermore, 
all these four scalars commute with (7 ■ 11)^; consequently, 

[M, (7 ■ n)2] = 0 . (3) 

It follows that the mass operator will be diagonal in the basis spanned by the eigenfunctions 
of (7 ■ Uf: 

- (7 • n)Vp(a;) = (4) 

If we work in the chiral representation in which S3 = i7i72 and 75 are both diagonal with 
eigenvalues a = ±1 and y = ±1, respectively, the eigenfunctions of (7 ■ fl)^ has the general 
form 


V’p(^) Fpcrx{.x)uJay-, (5) 

where are bispinors which are the eigenvectors of S3 and 75. The exact functional 
form of the Ep(j^{x) will depend on the specihc external held conhguration. Our method 
is based on the use of these eigenfunctions as basis functions. This is a natural choice as 
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they are the wavefunctions of the asymptotic states when a uniform external held is present. 
The advantage of using this representation is obvious: M can now be put in terms of its 
eigenvalues, so the problems arising from its dependence on the operator fl can be avoided. 

We now restrict our consideration to the case of a constant magnetic held of strength H, 
the vector potential of which may be taken to be = ( 0 , 0 , ifxi, 0 ), p = 0 , 1 , 2 , 3 . Our 
metric has the signature g^i, = diag(— 1 , 1 , 1 , 1 ). The eigenfunctions Ep^xi^) are now given 
by 

Ep^{x) = ( 6 ) 


where is a normalization factor and Dn{p) are the parabolic cylinder functions [Q with 
argument p = ^j 2 \eH\[xi — and indices (which are the quantum numbers of the Landau 
levels) 


n = n{k, a) = k + 


eHa 

2 \eH 


1 

2 ’ 


n = 0,1, 2,... 


( 7 ) 


Note that, in the absence of an external electric held, the eigenfunctions do not depend on 
y. The eigenvalue p stands for the set {po,P2,P3, k), where k is the quantum number of the 
quantized squared transverse momentum: 


- (7 ■ nx)Vp(a:) =-(7'ni + 7'n2)Vp(a:) 

= + eHT.^)^p{x) 

= 

= 2 \eH\kipp{x). (8) 

Notethat (7-n)2 = (7-n|,)2 + (7.n^)2, where (7-n|,)2 = (70no+73n3)2 = ng-H^, hence 
p^ = —pQ+pl + 2 \eH\k. The allowed values for k are seen from Eq.(^ to be fc = 0 , 1 , 2 , ... . 

Following Ritus 0 , we form the eigenfunction-matrices Ep{x) = diag(Fpii(a:), Ep_ii{x), 
Epi_i{x), Ep_i_i{x)). As noted above, in a pure magnetic held, references to y are irrelevant 
and can be dropped, hence 
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( 9 ) 


cr 

a 

Using the orthogonal property of the parabolic cylinder functions P, 

/ CO 

dpDn'{p)Dn{p) = V^nldnn’, 

-CO 

it is straightforward to establish that the Ep are orthonormal {Ep = 

j d^xEp/{x)Ep{x) = {27r)^S^^\p — p') 

= (27r)^4fc/(5(po - Po)5{p 2 - P2)S{P3 - Ps) 

as well as complete: 


( 10 ) 


( 11 ) 


'jld‘^pEp{x)Ep{y) = { 27 rY 6 ^^\x - y), j:d^p = Y, J dpodp2dp3 (12) 

k 

provided that the normalization constant in Eq.(P) is taken to be N{n) = (47r|eif|)^/‘^/\/rd. 
Since the Ep are linear combinations of the eigenfunctions of the mass operator, they satisfy 


d‘^x'M{x,x')Ep{x') = Ep{x)TiA{p)-, 


( 13 ) 


where S/i(p) represents the eigenvalue matrix of the mass operator. The Ep also satisfy the 
important property that 


7 ■ n Ep{x) = Ep{x) -f-p, 


( 14 ) 


where po = po, pi = 0 , p2 = —sgn(eif)Y^^ 2 |e^^, p^ = p^. Note that, due to the ro¬ 
tational symmetry about the direction of the magnetic held, the system is effectively a 
(2-1-1 )-dimensional one, as is evident from the momentum p. 

By using the above properties of the Ep-functions, the fermion Green’s function may be 
expressed as 


GA{x,y) 


d^p 


Ep{x) 


1 

7 - pT Sa(p) 


Ep{y), 


( 15 ) 
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Eqs.([T^) and (|T 3 |) guarantee that Eq.(|l|) is satisfied. It follows that, in the Ep-representation, 


Ga{.P,p') = J d'^xd‘^yEp{x)GA{,x,y)Ep>{y) 

= — p')- 


( 16 ) 


7 ■p + Sa(p) 

which shows explicitly that the fermion propagator is diagonal (in momentum) in this rep¬ 
resentation. Similarly, the mass operator may be written in the Ep-representation as 


M{p,p') = J d'^xd^x'Ep{x)M{x,x')Ep/{x') 

= (2lT'fS^‘\p - p')tA(p). 


(17) 


We may now write down the Schwinger-Dyson (SD) equation for the fermion self-energy. 
We shall work in the quenched, ladder approximation in which 


M[x,x') ^ ieA^GA{x,x')A- x'), 


( 18 ) 


where Dpy{x — x') is the bare photon propagator. 


(19) 


( 27 r)‘^ — ie 

The SD equation in this approximation reads 


{2n fS^‘^\p - p')Ea{p) d^xd^x'j: -^^Ep{x)YEi 


X 


Ep"{x)Y Epi{x)D pia{x - x). 


( 20 ) 


7-P"-FSa(p") ^ 

The SD Eq.(pOD can be simplified by performing the integrations over x, x', Pq, P2^ 

Pg exactly. Consider first the x-integrals. The Xq-, X2-, and Xa-integration each yields a 
^-function, leaving 

f d^xEp{x)YEp„{x)e^'i-^ = ( 27 r)=^( 5 ® (/ + g - p) ^ N{n)N{n”) 

(T, 

/ OO 

dxiD„(p)D„.(p")e*‘^i"i 7 °A 707 ^A", (21) 

-00 

where (p" + q-p) = 5 (po + qo -Po)8{p2 + ^2 - P2)8{p3 + gs -Ps), p” = ^j2\eH\{xi - ||), 
n" = n{k'\ cr"), and A" = A(ct"). Due to the presence of the 5(p2 + q2 — P2), the remaining 
xi-integral may be more conveniently written as 
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/ oo 1 


(22) 


where 


/ OC 

^^e*'?sgn(erf)gi^^(^ - q 2 )Dn"{r] + q 2 ), 

-CO 


(23) 


p = p + q 2 , and are dimensionless variables defined as 

q^,^J 2 \^\ 


g/2 = 


/i = 0, 1, 2, 3. 


2eH ’ 

If we transform to the polar coordinates: {q± = \Jqi + ql, ^ 
evaluate the Inn" by first noting that they satisfy @ 

dinn" {HA-i 9 ^) 


dp 


= isgn{eH){n - n)Inn"{qx, p)- 


(24) 

arctan(g 2 /gi)), we can 

(25) 


Hence, 


Inn"{q±,p) = 


where 


/ OO 

-CO 


Note that Inn"{,q±) = In"n{q±)- To compute Inn"{q±), we use the relation, 

mm{n,n") 


n\n"\ 


DMD„,.M = e-"'" , 

ml[n — m)l [n” — m )! 


—2m (h), 


and the Rodrigues formula. 


ri 

D4„) = 


to secure 


Inn" (gi) = ^2ire '"V„„"(gi), 


(26) 


(27) 


(28) 


(29) 


(30) 


where 
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min(n,n") 




nm 


ff\ 


m\{n — m)\{n''— m)\ 
Putting all the pieces together, we have 


[zsgn(eif)gj 


'\n-\-n''—2m 


( 31 ) 


1 




E 




( 32 ) 


\/n\n"\ 

where we have used the identity, 7°Ay*’ = A. Similarly, the x'-integrals yield 

I d^x'Ep.{x')YEpix')e-^'^-^' = ( 27 r)^ 5 ®(p//^^_p/)g-%K+p',)/( 2 erf) 

a',a" vnln I 

where n' = n{k',a'), n" = n{k",d"), A' = A((t'), and A" = A(a'"). The presence of the 
( 5 -functions in Eqs.(^) and (| 33 D allows easy integrations over p^, P2, and pg in Eq.(^), 
yielding 5 ^^\p — p') = 5 {po — Pq)6{p2 — P2)HP3 ~ Pa) which matches that on the left hand 
side of Eq.(l^). The SD equation is therefore reduced to 

/ pip. Asgn(eH){n-n"+h"-n')ip 


e-« Js"„-(9±)7 - (1 - 


Ay^A" 


-A"7^A', 


( 34 ) 


7 ■ p" Sa(p") 

where the summation over {(t} means summing over (T, a', a”, and d", and the momentum 
p" is understood to be: Pq = po — qo, p'l = 0 , P2 = — sgn{eH)^j 2 \eH\k", pg = pa — ga- 
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III. SOLUTION TO THE SCHWINGER-DYSON EQUATION 


An approximate solution to Eg. P^) may be obtained by observing that, due to the factor 
-2 

e“U in the integrand, contributions from large values of g± are suppressed. Thus, by keeping 
only the terms with the smallest power of g± in Jnn"(Q±), he., 

[max(n,n")]! ^in-n"i 

Jnn"(g±) —JZ —(*sgn(ehf)g±)l ' 

(35) 


\n — n" 
n\ 5r, 


and similarly for Jh"n'{,g±), the SD equation is simplihed to 

tA{p)6kk' ^ te\2\eH\)Y,Y.^^n"^n"n' f 

11 / r ^ yZilT j Q 

-A"7'^A'. 


k" {ct} 




AYA"- 


(36) 


7 ■ p" + Sa(p") 

The solution is expected to have the form S^(p) = Z{p)z ■ p + Ea(p), where Sa(p) is the 
dynamically generated fermion mass and is assumed to be proportional to the unit matrix 
(see remarks after Eq.(|^)). Eq.(^6D then reads 

e-9i 


[Z(p)7-p + S^(p)]4fc/ ~ Y{2\eH\)^^5nn"5n"n' j ■ 


where 


k" {.} (2^)^ 

T.a{P")[Gi - (1 - Qhhi] - [1 + Z{p")][G2 - (1 - OIE 2 ] 
[1 + Z{p')fp'^ + TPa{P') 


Gi = AYA"A'P,A, 

fUi ^ 4^(7-g)A"A"(7-g)A', 

qz 

G, = A7'^A"(7-/)A"7^A', 

IU 2 = ^A(7 ■ g)A"(7 ■ p")^''(7 ■ q )^'- 


(37) 


(38) 


The matrices Gi ,2 and lUi ,2 niay be simplified as follows. First we note that the A- 
matrices may be expressed as 




(39) 
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where = (^o-i — -D^ = 1. They also satisfy the commutation relations, 


A 7 I = 7^V 

and 


(40) 


[A.Tf] =0 = ^75*], (41) 

where the subscript T refers to the transverse components, p = 1,2, the subscript || refers 
to the longitudinal components, p = 0, 3, and V is the complement of A: 

V((r) = 1 - A(<t) = 1 (1 - C,Es). (42) 

Secondly, products of the A-matrices may be expressed in terms of a single A-matrix, e.g., 

A(a)A(a") = <5^a"A(a), (43) 


and similarly for the V-matrices. Using these relations, we hnd that, after performing the 
summation over the spin indices. 


^ ) ^nn'Aji"n'Gri ^kk" T ^k,k"—sgn{eH)A(^l.^ -|- 6 k^k''-i-sgn{eH)A(^ 1) 

^ ] ^nn'An"n'^^l — ^kkAkk") 

W} 

^ ) ^nn'An"n'^2 ^^kk' (h ' P-l) T T ’ (p|| 


W} 


sgn(eir)A(l) T ^k,k"+sgn{eH)^i, 1) i 


(44) 

(45) 

(46) 

^nn"Sn"n'W2 ~ 4fc'4fc" | 7 ' ^ (d ' ?||) 7 ' (P|| “ ?||) (7 ' ?||) [ • (47) 

{a} [ Q } 

where 7 ■ = 7 ^P 2 , A(l) = diag(l, 0, 1, 0), and A(—1) = diag(0, 1, 0, 1). In accordance 

with the small q± approximation used on the Jnn" and Jn>/n', Eg. (|35D , terms in lUi and IU 2 
that are proportional to q± have been dropped. It is satisfying that the spin summation 
produces the Kronecker delta, 6 kk', which matches the one on the left hand side of Eg. ( 0 )- 
Note also that, due to the restriction n = n" which arises from the small q±_ approximation, 
the summation over k” is now restricted to only three terms: for a given k, k" = k, k ±1. 
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The SD equation may now be written as 


z{ph -p + ^Aip) 


- - *e^(2|ei/|)^ 


e-^i 


(27r)4 g2 [1 + Z(p")]V + 

[1 + Z{p")] 5k"k ((1 + Ot ■ Pi. - (7 ■ ?||) 7 ■ {p\\ - ?||) (7 ■ ?||) 




T 2 T dk"^k—sgn{eH)^{ 1)^ 7 ' (p|| Q\\^ 

Z^A^P ) dk"k{^ T 1^) T 2 ^5fc",fc+sgn(er/)^(l) T sgn(eH)‘^( ^))] } ‘ 


(48) 


Recall that = {p\\ — gy)^ + 2\eH\k". 

Eq ■ (|48|) shows that our earlier assumption of Eyi(p) being proportional to the unit matrix 
is correct only for the k" = k term. The reason is that 7 ■ If, and hence the mass operator 
M, does not commute with S 3 . However, if we consider only the low energy {p^ -C |ei7|) 
behaviors, in particular, in the p± = Q = k limit, the k" = 0 term will dominate and we 
obtain the approximate SD equation. 


Z{p\\)l-P\\ +^a{p\\) 

~ — ie^(2|eTr|) j 


Q-ql 


1 + Z(p|| -gy) 


(2vr)4 g2 [1 + Z{p\\ 
1^-1 


Til -9||)^ + Si(py -gy) 

(7 ■ ^ 11 ) 7 ■ (p|| - gy) (7 ■ ?||) + ^a{p\\ - gy)(l + 0 f • (49) 


In this case Ha is proportional to the unit matrix. This approximation is equivalent to the 
lowest Landau level approximation employed in |^. 

We see from Eq.(^UD that, in the Feynman gauge (.^ = 1), ^(py) = 0 and the SD equation 
for the dynamically generated fermion mass becomes 

d^g e-«i SA(py-gy) 


HAiPii) ^ e\A\eH\) 


(50) 


(27r)4 g2 (pii _ ^11 )2 + _ ^ 11 ) ’ 

where we have made a Wick rotation to Euclidean space: po —>• ip^, go —> iq^. This can 
be turned into a differential equation for Ha, as was done in Ref. where the momentum 
dependence of the dynamical mass is discussed. Here we content ourselves with Ending a 
solution for the infrared fermion mass scale. 
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Ea(0) ~ e^{MeH\) 


d‘^q e ^^(gii) 


(51) 


(27r)4 g2 2|eff|gJ + S^(g||)- 
Since Sa(q'||) is expected to diminish with increasing the integral in Eq. (|5l|) is dominated 
by the contribntions from small gy. It is therefore reasonable to approximate S^(g||) in the 
integrand by S^(0) = m x 1, where m is the dynamical mass and 1 is the nnit matrix, thns 
securing the gap equation, 

e-^i 1 


1 2/.I TTW f “ y 6 

“ ® ^ J (27r)4 g2 2\eH\q^ + m^ 


47r^ 


leifl 


1 


ffl 

poo poo 0 

L Jo + gy 2|ehf|gj| + 


a I ln(2|eif|g^/m^) 

^ 2|ei7|gi-m2 ' 


The solution to Eq.(|5^ has the form 


(52) 


m ~ a ^/\^\ 


(53) 


where a and b are positive constants of order 1. The nonperturbative nature of this result 
is apparent. Furthermore, according to the last equality of Eq.(^), the dominant contribu¬ 
tions to the integral come from the region 2\eH\qj_ ~ Consistency with our small q± 
assumption requires that m -C ^J\eH\, which in turn requires that a 1. In other words, 
the solution for the dynamical mass found above applies to the weak-coupling regime of 


QED. 

To establish that the above solution to the SD equation for the fermion self-energy does 
indeed correspond to a dynamical chiral symmetry breaking solution, it is necessary to 
demonstrate the existence of the corresponding Nambu-Goldstone (NG) boson. One way to 
establish this is by studying the Bethe-Salpeter equation of the bound-state NG boson, as 
was done by Gusynin et al. [^, who found a solution consistent with our Eg. (|53|) . We show 
in the Appendix how the same solution, Eq. (|5^) , can be obtained from the Bethe-Salpeter 
equation for the NG boson, using the Ep-representation of the fermion progagator. This 
helps to justify that an external magnetic held serves as a catalyst for chiral symmetry 
breaking in QED. 


11 










IV. SYMMETRY BREAKING AT NONZERO TEMPERATURE 


The formalism described above is very useful for studying dynamical symmetry breaking 
in an external field at nonzero temperatures. As expected on physical grounds 0 , the long- 
range order of a system decreases as the temperature increases and chiral symmetry will 
generally be restored above a certain critical temperature. We shall see that this expectation 
is indeed correct for the case of chiral symmetry breaking in a magnetic held and an estimate 
of the critical temperature will be obtained. 

To incorporate the thermal effects, we use the imaginary time and energy formalism 
in which 


0 < ixQ < (3 


(54) 


and 


qo = 21 'tt for bosons, /'= 0, ±1, ±2,... 

Po = (2/+ l)vr/(—i/3) for fermions, / = 0, ±1, ±2,... (55) 


where (3 = 1/T, with the Boltzmann constant ks = 1. The analysis given in Sections II 
and III can be repeated for the case of hnite temperature by implementing the following 
replacements 0: 


d3^x 

4^ 


r-i/3 


:d^p 

— p') 


Jo 

27T 

—1(3 


dxo j d^x 


/ dp2 / dp3 


(56) 


l k 


S?{p -p) = -i^dw5kk'd{p2 - P2)^{po - P3) 


-i(3 


271 


and for the photon propagator 


d'^q 


The gap equation now reads 


27r 


—i(3 


^J2 dqi dq2 / dq^ 


(57) 


2 (y 

1 ~ —T\eH\ 

7T 


('OO POO 

dqs / dqle~ 

—00 JO 


E 


rQ 2 + Qi + 7r2T2(2/' - 1)^ 


(58) 
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where Qi = ql + rn^ and Q 2 = q^ + 2\eH\qj_. Following Ref. [jl^, we snm the series by means 
of the Poisson snm formnla, which states that, if c(r) is the Fonrier transform of 6(0;), 


/ CO 

-CO 


the following identity will hold: 


Here, taking 


E KH = E c(27rA). 

'=—00 A=— CO 


(69) 


(60) 


b(ui) = 


Q 2 + 4 HW- Qi + 7 iV^{2u: - 1)'^’ 


(61) 


we have 


c r = 


87rT2 J-c 


duexp 


u \u\ 

+ 




(62) 


2 2nT'‘ 2nT 

After integrating over n, we make nse of the Poisson snm formnla to rewrite the snmma- 
tion in Eq.(^) in terms of a more manageable series: 

1 1 


E 


Q 2 + 47r2T2/'2 + 7r2T2(2/' - 1)2 

1 1 


2Tv/^,5 l 


e I’ + (—l)^e T 


\fQ\ + \fQ2 


4 {"^Qi + V + vr2T2 


e T _ (-l)-^e T ) 




(63) 


'' iVOi ~ _ 

where the prime on the snmmation sign means that the A = 0 term is connted with half 
weight. Snbsitnting Eq. (|6^) into Eq. (|^) and snmming over the inhnite geometric series, we 
obtain 


01 . . i 

1 ~ —|ehf| / 

TT J —00 


roo 

dqs / 

Jo y/QlQ 2 


+ 


y/Qi + VQ 2 


_VS 2 . -1 4 /^ 

.1 — e T 1 + e T 


{y/Qi + a/ Q 2 Y + 
y/Qi — 


_V^ 


1 — e T 1 + e T 


(y^_y^)2 + 7r2T2 


(64) 
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Note that the only approximations we have made so far are the quenched, ladder approxima¬ 
tion (Eq.(^)), the small q± approximation, and keeping only the k" = 0 terms in Eq.(^). 
Aside from these approximations, the hnite-temperature gap equation, Eq.(^), is exact in 
its dependence on the coupling constant, the magnetic held, and the temperature. 

We shall consider below the zero, low, and high temperature limits of Eq. (p^) . For T = 0, 
Eq.(|64|) is reduced to 

1 


/y POO POO 

1 ~ -\eH\ / dqs / dq 

TT J — OO J 0 


~ 2 


(65) 


\fQ\Q 2 \fQ\ + \fQ2 

which, one can easily check, is just what Eq. (|5^ becomes when the integration there over 
g4 is done. Thus, we have recovered the T = 0 result. 

For the low temperature case, we can approximate the sum over /' in Eq. (|^ by the A = 0 
term (with half weight) on the right hand side of Eq. (|63D , the other terms being exponentially 
small. Treating the thermal effects as a perturbation, we write m?{T) = -|- 5m? with 

5m? « ttIq, where ttiq is the fermion dynamical mass for T = 0 (i.e., the solution in 
Eq.(^). At small T, the difference between Eq.(|5BD and Eq.(|B3D yields 

( 66 ) 




J_2 

h 


where 


/ OO POO 

dqj dgie-^x 

-OO ^0 

/ OO POO 

dq, / dqle-^^ 

-OO 


1 1 


+ 


h = 


Q \/Q 2 \fQ + V Q2 \ \fQ y/Q + VQ2 / 

1 1 1 


(67) 


v/Q (TQ + 

with Q = rriQ -|- g|. Note that the infrared regions of q, and qj_ dominate the integrals, 
just as in the T = 0 case. More importantly, both Ji and I 2 are hnite and positive so that 
5m‘^ is negative (and small, being proportional to T^). Thus, thermal effects tend to reduce 
the fermion dynamical mass, i.e., chiral symmetry tends to be restored as the temperature 
increases, as discussed earlier. However, we should stress that, so long as T is small, a 
non-zero dynamical mass will be generated in the presence of the magnetic held. 
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At high temperatures (T > ^J\eH\), the /' = 0 term dominates the sum in Eq.(^8|), thus 
yielding 


1 ~ — T\eH\ / dgs / dq‘^e~^^^- 

7T J —CO J 0 


( 68 ) 


Q2Qi + 7T^T^' 

Since the dominant contributions to the ga-integral come from small values of g|, we may 
approximate the denominator (Qi + tt^T^) by [m^(T)+7r^T^]. After evaluating the integrals, 
we obtain 


rn?{T) + tt^T^ ~ aT\J2n\eH\ (69) 

Obviously, at weak couplings (a -C 1) and at high temperatures (T > ^J\eH\), there is no 
non-negative solution for m?{T), i.e., no chiral symmetry breaking solution. Thus, both the 
high temperature and the low temperature solutions indicate that chiral symmetry will be 
restored at high temperatures. 

If we define the critical temperature for chiral symmetry breaking to be the temperature 
at which m^{T) vanishes, i.e., m‘^{Tc) = 0, Eq.(|69[) provides an estimate of the critical 
temperature: 


a 


- ^v27r|eif| 


TT^ 


~ 1.5 X 10“2 eV 


\H\ 


(70) 


gauss 

For a magnetic field of any given strength, the critical temperature can be more exactly 
calculated by numerically solving the gap equation, Eq. ([0^) , with Qi replaced by g|. 

Finally, we consider the suggestion by Gusynin et al. that the chiral symmetry break¬ 
ing solution found above may play a role during the electroweak phase transition in the early 
universe. The electroweak phase transition took place at a temperature of order 100 GeV. 
From Eq.(^), this requires a magnetic field of 10^° gauss or stronger, which is significantly 
larger than any estimates of the magnetic field strength at the time of the electroweak phase 
transition El- We therefore conclude that the chiral symmetry breaking solution considered 
here does not play any role in the electroweak phase transition. 
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V. CONCLUSION 


We have described a formalism for studying the physics of chiral symmetry breaking in 
an external held via the nonperturbative Schwinger-Dyson equation. The T^p-representation 
for the fermion propagator proposed here has the advantage that the dependence on the 
operator fl^ is removed. It also has the advantage over Schwinger’s proper time formalism 
for calculating hnite-temperature effects. 

We have applied our method to examine chiral symmetry breaking in QED in a constant 
external magnetic held. We hnd that, even when the coupling constant is small, an external 
magnetic held can trigger the dynamical breaking of chiral symmetry in QED, with the 
dynamical mass of the fermion given by Eq. (|^ . Our result agrees with that found in 
Ref. @] which used an approach rather diherent from ours. The existence of the Nambu- 
Goldstone boson of chiral symmetry breaking is demonstrated in the Appendix by showing 
that the same solution (Eq. (|53|)) solves the Bethe-Salpeter equation for the bound-state 
Nambu-Goldstone boson. 

We have also obtained an estimate of the critical temperature Tc of the aforementioned 
chiral symmetry breaking, Eq.(^). Ghiral symmetry is a good symmetry above this critical 
temperature, and it is spontaneously broken at temperatures below T^. Our result ren¬ 
ders invalid the suggestion in the literature that the chiral symmetry breaking solution 
described above may be relevant for the electroweak phase transition in the early universe. 

There remain several interesting questions which we are investigating. For instance, how 
do other background field configurations affect chiral symmetry breaking? (As an example, 
it is expected that an electric held would tend to break up the condensate and destabilize the 
vacuum, thus inhibiting chiral symmetry breaking |T^.) Are there strong-coupling solutions 
of chiral symmetry breaking in an external magnetic held? What are the ehects of additional 
four-fermion interactions? Recall that four-fermion operators play a crucial role in obtaining 
a consistent chiral symmetry breaking solution in quenched, ladder QED [|T^. We hope to 
report our hndings in the near future. 
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APPENDIX A: NAMBU-GOLDSTONE BOSONS 


We verify in this appendix that the dynamical mass found in Eg .([531) corresponds to a 
spontaneous chiral symmetry breaking solution by examining the Bethe-Salpeter equation 
of the Nambu-Goldstone (NG) boson. In the quenched, ladder approximation, the Bethe- 
Salpeter equation of the NG boson has the form [Q 

0(x, y;P) = - ie^ J dV J d^y'GAix, x')'^^4>{x', y'; P)YGA{y', y)Df,^{y’ - x'). (Al) 

In the Ep-representation, this can be expressed as 

(27ryS^^^(p-p')(^(p;P) = - ie‘^ ( d^xd^yff^^ -^— Ep{x)YEp>'{x) 

J J 7 ■ p -|- m 

• 0(/; P)Ep'iy)YEp'{y) —— Dp^{y - x), (A2) 

7 ■ p -|- m 

where 

0(x, y;P)=p -^^Ep{x)^{p; P)Ep{y) (A3) 

and m is the dynamically generated fermion mass. 

If we introduce 


X(p; P) = (7 ■ p -h m) 0 (p; P) (7 ■ p m) , 


(A 4 ) 


we hnd 


{2tiY5^^\p - p')x(p; P) = - ie^ J ^^Ep{x)j>^Ep.{x) 


(2vr) 

-x{p";P)- 


^ ■ p" + m 7 ■ p" + m 

■ Ep"{.y)YEp'{.y)Dpu{.y - x). 

Note the similarity of this equation to Eg. (PU]) . After integrating over x and y (see Eq.i 
and Eq.(^)) as well as over Pq, P2 and P3, we obtain (in the Feynman gauge) 


(AS) 


x{p-,P)3kk’ = - ie^(2|eP|)^^ 


k" {<t} ' 


(2w) 


vrd 


n'\n"\n"\ 


Q-ql 


Jnn" iGx)Jn”n' (?±) A 


1 


■y ■ p" + m 


x{/;P) 


X ■ p" + m 


A"7mA', 


(A6) 
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where = po + go, K = 0, P 2 = - sgn{eH)^j2\eH\k", pg = ps + go- 

We are interested in the P = 0 behavior of x{p] P)i which is expected to have the form 


X(p; 0) = A(p)75, 

where A{p) is a scalar function. It follows from Eg. (|A6|) that 

r p-*sgn(eH)(n-n"+n"-n')(/j 

A{p)5kk’ = ie^(2|ehf|) ' 


(A7) 


e-^i 


k" {a}'' Vn\n'\n”\n”\ 

Jnn"[(l±)Jh"n'[Q±)' 


- - - - - ' p>/2 _|_ ^2 ’ 

with Gi dehned in Eg. (|38|) . After taking the small q± approximation, Eq.( 
simplihed to 

A{p)6kk' — ( 2 1 eif I ) ^nn" ^h"n' / /„ \4 ^2~W2~T 2‘ 

■> (27r)4 g^ p''^ + m^ 

If we consider the infrared behavior of A(p), we have 

d^q e'-^i A(g||) 


(A8) 
this is 

(A9) 


A(0) ~ — ie^(4|eih|) J 


(AlO) 


(27r)‘^ g2 2|ei7|gj| + ’ 

where we have performed the spin summation and ignored the k" = 1 term, in accordance 
with the approximation used in obtaining Eg. (|^) . Transforming to Euclidean space and 
noting that the integral is dominated by contributions from small gy so that A(g||) in the 
integrand can be approximated by A(0), we secure 

d'^q e“'^J- 1 


1 ~ e^(4|ei7|) J ■ 


(27r)^ g2 2|ei7|g|| + ’ 

which is the same as the gap equation. Eg. (|5^ . Thus, the dynamical mass found in Eg. 
provides a consistent solution to the Bethe-Salpeter equation for the NG boson. 


(All) 
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